Proximity induced interface bound states in superconductor-graphene junctions 
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We show that interface bound states are formed at isolated graphene-superconductor junctions. 
These states arise due to the interplay of virtual Andreev and normal reflections taking place at 
these interfaces. Simple analytical expressions for their dispersion are obtained considering interfaces 
formed along armchair or zig-zag edges. It is shown that the states are sensitive to a supercurrent 
flowing on the superconducting electrode. The states provide long range superconducting correla- 
tions on the graphene layer which may be exploited for the detection of crossed Andreev processes. 

PACS numbers: 73.23.-b, 74.45. +c, 74.78.Na, 73.20.-r 



Introduction: Several striking transport properties 
have been predicted to emerge from the peculiar elec- 
tronic band structure of single atom graphite layers, 
known as graphene [l| . Of particular interest is the case 
of a graphene layer in contact with a superconducting 
electrode, a situation which has been explored in several 
recent experiments @|. Here, as in the case of normal 
metals, the mechanism which dominates the electronic 
transport at subgap energies is the Andreev reflection, 
i.e. the conversion of incident electrons into reflected 
holes with the creation of Cooper pairs in the supercon- 
ductor. However, while in the case of normal metals the 
reflected hole has typically the opposite mean velocity 
to the incident electron (retroreflection) in the case of 
graphene it is possible to have Andreev reflections with 
a specular character, as first shown in 

Transport and electronic properties at graphene- 
superconductor junctions have been analyzed in several 
works [i[ . ft has been shown that the special character of 
Andreev reflection in graphene leads to modifications in 
the differential conductance compared to that of conven- 
tional N-S junctions [1,0]. The effect on the local density 
of states (LDOS) has been studied in Refs. d, 0- The 
properties of bound states arising from multiple Andreev 
reflections in S-graphene-S junctions have been analyzed 
in Refs. Q. However, as we show in this work, the spe- 
cial electronic properties of graphene are such that bound 
states can be formed even at isolated single junctions. 

The mechanism for the emergence of these states can 
be understood from the scheme depicted in the left panel 
of Fig. [TJ As is usually assumed the junction can be 
modeled as an abrupt discontinuity between two regions 
described by the Bogoliubov-de Gcnnes-Dirac equation, 
taking a finite superconducting order parameter A and 
large doping Ep > A on the superconducting side and 
zero order parameter and small doping Ep ~ A on the 
normal side. For the analysis it is instructive to include 
an artificial intermediate normal region with A = and 
Ep = Ep, whose width, d, can be taken to zero at the 
end of the calculation. This intermediate region allows 
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FIG. 1: (Color online) Simple model for the emergence of 
IBSs (left panel). It illustrates the scattering processes tak- 
ing place at a graphene-superconductor interface with an in- 
termediate heavily doped normal graphene region of width d. 
Cases (i) and (ii) correspond to the case hvq < \E ± Ep\ and 
hvq > \E±Ep\ respectively with A > E > Ef- (Right panel) 
Graphene-superconductor junctions along different edges. On 
the superconducting side (shaded areas) the on-site order pa- 
rameter A is finite and the doping level is high (Ep S> A). 



to spatially separate normal reflection due to the Fermi 
energy mismatch from the Andreev reflection associated 
to the jump in A. As shown in Fig. [I] (case i), an incident 
electron from the normal side with energy E and parallel 
momentum hq such that hvq < \E— Ep\ is partially trans- 
mitted into the intermediate region and after a sequence 
of normal and Andreev reflections would be reflected as a 
hole. This process can either correspond to retro or spec- 
ular Andreev reflection depending on whether E < Ep 
or E > E F Q. For hvq > \E ± Ep\ neither electron 
or holes can propagate within the graphene normal re- 
gion. However, virtual processes like the one depicted in 
Fig. Q] (case ii) would be present. These correspond to 
sequences of Andreev and normal reflections within the 
intermediate region. A bound state emerges when the 
total phase 4> accumulated in such processes reach the 
resonance condition d> = 2nir. 
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The aim of the present work is to demonstrate the 
existence of these interface bound states (IBS) and to 
analyze their properties for different types of graphcnc- 
superconductor junctions. After completing the anal- 
ysis for the case of the simple model sketched above, 
which implicitly assumes a decoupling of the two val- 
leys in the graphene band structure, we consider more 
microscopic models for junctions formed along armchair 
or zig-zag edges. We study the effect of an additional 
potential barrier at the interface and the possibility to 
modify the states by a supercurrent flowing through the 
superconductor. We finally discuss the potential use of 
these states for the generation of non-locally entangled 
Andrccv pairs. 

It is quite straightforward to determine the dispersion 
relation for the IBS from the model represented in the left 
panel of Fig. [TJ The phase accumulated by a sequence 
of normal and Andreev reflections in the intermediate re- 
gion can be obtained from the corresponding coefficients 
r e , rh and ta- Following Ref. Q one obtains 
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arcsm hvq/(E±E < p>). The condition 
Ep ^> A, E, hvq allows to take ot\ h ~ 0. On the other 
hand, in the region of evanescent electron and hole states 
for graphene (|ft«(7| > \E±Ep\) r e ^ become a pure phase 
factor e %tpB > h , withi/?^ = — 2sigs\(q/ (E±Ep)) arctane Ac -' 1 
and A e ,h = sign(<j)arcosh(ft«(7/|i? ± Ep\). For the An- 
dreev reflection coefficient between regions / and S one 
has ta = & lVA where ipA = arccosS/A, as it corresponds 
to the Andreev reflection at an ideal N-S interface with 
Ep ^ A [i| . In the limit d — > the total phase accumu- 
lated is thus <p = 2<PA + <Pe + <Ph, from which one obtains 
the following dispersion relation 



E e (A c +A h )/2 



A 



= ±- 



sign(i? 2 



2-^/cosh A e coshA/j 



(2) 



This dispersion simplifies 



to 



E/A 



zthvq/ (hvq) 2 + A 2 at the charge neutrality point 
(i.e. for Ep = 0). In this case the IBS approaches 
zero energy for q — > and tend asymptotically to 
the superconducting gap for large q. Notice also that 
the decay of the states into the graphene bulk region 
(x < in the left panel of Fig. [TJ is set by e x ^"- h , 
where £ e ,fc = f>-v/(\E ± Ep\smh(X e ^)) for the electron 
and hole components respectively, which can be clearly 
much larger than the superconducting coherence length 
£o = frv/A when Ep <C A. It is also interesting to notice 
that the IBSs survive when Ep > A, i.e. in the regime 
corresponding to the usual Andreev retroreflection, but 
with a much smaller spatial extension. 



In order to analyze the existence and the charac- 
teristics of the IBSs for different types of graphenc- 
superconductor junctions we make use of the Green func- 
tion formalism based on tight-binding models for these 
junctions which was introduced in Ref. Within this 
formalism the retarded green functions at the interface 



G(E, q) are given by 



where g corresponds 
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to the surface of the uncoupled semi-infinite graphene 
layer and S is the self-energy associated to the coupling 
with the superconductor. In general all these quantities 
have a 2 x 2 structure both in the sublattice (indicated 
by the hat symbol) and the Nambu (indicated by the 
check symbol) spaces. Once these quantities for each 
type of interface have been determined, the existence 
of an IBS can be established by analyzing the equation 
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Interface along an armchair edge: We first consider 
an interface constructed along an armchair edge, as 
schematically depicted in the right panel of Fig. [TJ In 
a rather generic way one can write g = g e {fo + f z )/2 + 

9h{fo - f z )/2 and £/i 9 = Pf z g B csf z + jf z a x , where 
g e ,h describe the propagation of e and h components in 
the uncoupled graphene layer, gscs = 9% + f^x with 
g = -Ef /A = -E/VA 2 - E 2 being the BCS dimen- 
sionless Green functions, and (3 and 7 are parameters 
which allow to control the transparency and the type of 
interface. As discussed in Ref. 7 = corresponds 
to a model in which the coherence between the sublat- 
tices of graphene is broken on the superconducting side 
(bulk-BCS model), whereas for (i = v3/2 and 7=1/2 
corresponds to the ideal case where superconductivity is 
induced on the graphene layer by a superconducting elec- 
trode deposited on top, thus leading to a heavily doped 
graphene superconductor (HDSC). 

To make further analytical progress we take the limit 
E,A,hvq <C t g in g e h of Ref. Q, where t g denotes 
the hopping element between neighboring sites in the 
graphene layer. In this case and for hvq > \E± Ep\, g e ^ 
adopt the form t g g e< h = -\ (/i e ,fcOo + Ve,h&y) ± &x\ , 
where \i e _] % = sign(<7)/sinhA e ./ l and v e _] % = sign(E ± 
Ep )/ tanhA e ,/j. The Green functions matrix has the 



property g e \ 



-tigl \ =F t g (T x . Using this property and 



the definition for the self-energy the equation for the IBSs 
in this case becomes 



det [t 2 g gh9e + Pgt g {g e 
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For the HDSC model (i.e. /3 = \/i/2 and 7 = 1/2) the 
equation for the IBSs reduce to the one already found 
within the simple analytical model (Eq. j2])). This leads 
to a single root for arbitrary doping which is four-fold de- 
generate due to valley and spin symmetry. Fig. [5] shows 
a color-scale plot of the spectral density at a distance 



FIG. 2: Gray-scale plot of the spectral density at a distance 
~ £o from the interface defined along an armchair edge. The 
results were obtained using the HDSC model of Ref. for 
E F = (left panel) and E F = A/2 (right panel). The full 
fines indicate the position of the IBS determined from Eq. 

©■ 



~ £o from the interface on the graphene layer with two 
different doping conditions. The full lines correspond to 
the IBS dispersion obtained by solving Eq. @. As can 
be observed, the minimal energy for the IBSs, E rn i n , de- 
pends on Ep. Further analysis of Eq. ((5J) reveals that it 
satisfies the cubic equation E^ nin + E^^Ep — A 2 E F = 0, 
thus evolving between and A as Ep increases. The 
transition between E m i n > Ep and E m i n < Ep occurs 
at Ep = A/\/2. The presence of the IBSs manifests also 
in the appearance of singularities in the LDOS around 
E = ±A (see Ref. Q). The behavior of the LDOS is 
analyzed in more detail below. 

On the other hand, for the bulk-BCS model (i.e. 7 = 
and (3 G (0, 1)) one obtains 

-/?VOe/M + VeVh - 1) + V3f3g(n e , + + 
+ ^(l + 2^) + 3 (1 + 20 > ){UeVh _ ^ h) + 1 = . (4) 

In this case the degeneracy associated to the two val- 
leys in the band-structure of graphene is generally broken 
(except for Ep = 0). The roots gradually evolve towards 
the linear dispersion \E + Ep\ = hvq as [3 — > 0, which 
corresponds to the armchair edge state of the isolated 
graphene layer [Tol ] . 

Interface along a zig-zag edge: We now consider an 
interface along a zig-zag edge as illustrated in the right 
panel of Fig. [TJ The Green functions for the semi-infinite 
zig-zag edge can be obtained following the same formal- 
ism as in Ref. Q- In the continuous limit <? e ,/i becomes 



le ia e ,h zp e 
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where as in Eq. |T]) sin a e ,/i = (hvq)/(E±Ep) but with q 
measured with respect to the point K = 27r/3a, where a 
is the lattice constant indicated in the right panel of Fig. 
[TJ There exists an additional branch where q is measured 
from the opposite Dirac point at — K. The self-energy 



FIG. 3: (Color online) Dispersion relation for the IBSs on a 
zig-zag interface for decreasing parameter /3 controlling the 
coupling with the superconductor. The parallel momentum q 
in Eq. ([6]) is measured from the Dirac points at K = ±27r/3a. 



due to the coupling with the superconductor is in this 

case £ = f3t g (ao + o z ) {f z gBCST z ) /2. The equation for 
the IBSs then becomes 



E _ ± e (A e +A h )/2 _ sign ( g 2 _ E 2 p )/3 2 e -(A e +A h )/2 
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which looks very similar to Eq. ^ except for the presence 
of the parameter (3 controlling the coupling and the al- 
ready mentioned redefinition of the parallel momentum q. 
An interesting property of zig-zag edges is the presence of 
zero energy states for total parallel momentum between 
(— K, K) and Ep = When the coupling to the su- 
perconductor is turned on by increasing the parameter f3, 
one observes that the zero energy states evolve acquiring 
a finite slope. These states can thus be identified with the 
IBS for this type of interface. This is illustrated in Fig. 
[3] When the coupling parameter (3 reaches 1 the usual 
dispersion of the simplest analytical model is recovered. 

Effect of a supercurrent: a supcrcurrent flowing on the 
superconducting side of the junction modifies the spa- 
tial variation of the phase of the order parameter which 
produces a Doppler shift in the energy of the quasi- 
particles. This shift, obtained from the Bogoliubov-de 
Gennes-Dirac equations within the Andreev approxima- 
tion, is given by r\ = (hv) 2 q s q/Ep, where Hq s is the mo- 
mentum of the Cooper pairs assumed to be parallel to 
the interface. This result is equivalent to the one found 
in Refs. [12( for conventional and two-band superconduc- 



tors. Notice that for this analysis we go beyond the limit 
Ep — > oo taken in the initial simple model. The expres- 
sion for the reflection coefficients of Eq. (flj still holds 
but ckg ~ — a ! h = a 1 = axesin fivq/Ep is kept finite. On 
the other hand, the phase of the Andreev reflection coef- 
ficient between the intermediate region and the current- 
carrying superconductor becomes tpA = arccos£"/A(g s ), 
where E' = E + r\. At zero temperature, due to Lan- 
dau criterion, the order parameter is unaffected by the 
supercurrent while Hvq s < A (0) [12]. Therefore, in this 
condition, A(q s ) ~ A(0) = A. For the case E > Ep 
we thus get the following modified equation for the IBSs 
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FIG. 4: (Color online) Effect of a supercurrent flowing on 
the superconducting electrode on the dispersion relation (left 
panel) and on the local density of states at a distance ~ Co/ 10 
from the interface normalized to po = A(a/hv) 2 /2n (right 
panel). The results correspond to hvq s /A = 0.0, 0.25, 0.5 and 
0.75 with E% = 100A. 

within the simple model sketched in Fig. Q] 

E_ _ ± sinh (A e + X h )/2 + sin a 1 sinh (A e - X h )/2 
A yj (cosh A e — sin a 7 ) (cosh Xh + sin a 7 ) 

Figure |4] illustrates the effect of a supercurrent both in 
the dispersion relation of the IBS (left panel) and in the 
local density of states close to the interface (right panel). 
For q s = 0, the IBS manifest in a finite LDOS for E < A 
and a sharp peak at E = A. Qualitatively, the presence 
of a supercurrent breaks the symmetry with respect to 
inversion of the parallel momentum hq and leads to a 
splitting of the singularity at E ~ A in the LDOS. Note 
that this implies the appearance of an induced net cur- 
rent on the graphene side (for \x\ < £). For Ep = the 
distortion of the dispersion relation for finite and small 
q s is given by E(q s ,q) = E(0, q) + (hvq) 2 rj/ {(hvq) 2 + A 2 ). 

A more quantitative analysis of the effect of a super- 
current requires the estimation of the parameter Ep. 
This parameter is very much dependent on the fabri- 
cation methods and material properties of the metallic 
electrodes deposited on top of the graphene layer. Ac- 
cording to the ab-initio calculations of Ref. [l3| for Pd on 
graphene a typical estimate would be Ep ~ Q.leV, which 
for a superconductor like Nb gives a ratio Ep/ A ~ 100. 
The results on Fig. [4] have been obtained for this ratio. 

Conclusions: We have shown that interface bound 
states appear at graphene-superconductor junctions. 
The properties of these states are sensitive to the type 
of edge forming the interface, its transparency and the 
doping conditions of the graphene layer. We have demon- 
strated that the interface states evolve towards the edge 
states of the isolated graphene layer when the trans- 
parency of the interface is reduced. We have also shown 
that they can be modulated by a supercurrent flowing 
through the superconductor in the direction parallel to 
the interface. Even when our analysis has been restricted 
to interfaces along armchair or zig-zag edges we expect 
the appearance of IBSs to be a general property of any 
edge orientation. We also notice that inclusion of weak 



disorder along the interface introducing a small uncer- 
tainty in the parallel momentum dq would not prevent 
the emergence of IBSs provided that 5q <C A/hv. 

As a final remark we would like to comment that the 
existence of IBSs induce long range superconducting cor- 
relations between distant points on the graphene layer 
that are close to the interface. This property could be 
exploited to detect crossed Andreev processes and there- 
fore entangled electron pairs using weakly coupled STM 
probes on a graphene-superconductor junction, in a con- 
figuration like the one proposed in Ref. [3] ■ The analysis 
of non-local transport in this system will be the object 
of a separate work. 
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